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Abstract 

Background: Nuna bean is a type of ancient common bean {Phaseolus vulgaris L) native to the Andean region of 
South America, whose seeds possess the unusual property of popping. The nutritional features of popped seeds 
make them a healthy low fat and high protein snack. However, flowering of nuna bean only takes place under 
short-day photoperiod conditions, which means a difficulty to extend production to areas where such conditions 
do not prevail. Therefore, breeding programs of adaptation traits will facilitate the diversification of the bean crops 
and the development of new varieties with enhanced healthy properties. Although the popping trait has been 
profusely studied in maize (popcorn), little is known about the biology and genetic basis of the popping ability in 
common bean. To obtain insights into the genetics of popping ability related traits of nuna bean, a comprehensive 
quantitative trait loci (QTL) analysis was performed to detect single-locus and epistatic QTLs responsible for the 
phenotypic variance observed in these traits. 

Results: A mapping population of 185 recombinant inbred lines (RILs) derived from a cross between two Andean 
common bean genotypes was evaluated for three popping related traits, popping dimension index (PDI), expansion 
coefficient (EC), and percentage of unpopped seeds (PUS), in five different environmental conditions. The genetic 
map constructed included 193 loci across 12 linkage groups (LGs), covering a genetic distance of 822.1 cM, with 
an average of 4.3 cM per marker. Individual and multi-environment QTL analyses detected a total of nineteen 
single-locus QTLs, highlighting among them the co-localized QTLs for the three popping ability traits placed on 
LGs 3, 5, 6, and 7, which together explained 24.9, 14.5, and 25.3% of the phenotypic variance for PDI, EC, and 
PUS, respectively. Interestingly, epistatic interactions among QTLs have been detected, which could have a key role 
in the genetic control of popping. 

Conclusions: The QTLs here reported constitute useful tools for marker assisted selection breeding programs 
aimed at improving nuna bean cultivars, as well as for extending our knowledge of the genetic determinants and 
genotype x environment interaction involved in the popping ability traits of this bean crop. 
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Background 

Popbean or nuna bean (Phaseolus vulgaris L., Fabaceae) 
is traditionally grown in the Andean highlands of South 
America at 2,000-3,000 meters above sea level, where it 
is commonly sold in local markets or consumed at 
home, and it is thought to be an ancient pre-ceramic 
landrace [1]. It seems probable that nuna beans origi- 
nated in the Andes, where they are sympatric with wild 
and primitive common bean populations in certain parts 
of Peru and Bolivia, and they may have been present in 
the early stages of Andean agriculture [2,3]. The first se- 
lection pressures leading to domestication of common 
bean could have resulted in the development of popping 
beans, and it appears that toasting grains was a well- 
established tradition in the Andes and possibly in 
Mesoamerica, where early maize races have also been 
used for popping. However, no evidence of nuna beans 
has been found in Mesoamerica, most likely due to gen- 
etic differences between the Mesoamerican and Andean 
gene pools [1]. This, among other factors, may explain 
their contrasting popping ability and photoperiod re- 
sponse [4]. 

Popbean is tropical in appearance, with aggressive in- 
determinate type 4 growth habit [5] and day-length sen- 
sitivity since it requires nights of at least 11-12 h for 
flowering induction. Growing nuna bean in temperate 
areas requires the development of cultivars that are 
insensitive to photoperiod, like the modern dry bean 
cultivars. Furthermore, with a view to commercial pro- 
duction, determinate growth habit is desirable, since it 
reduces growing and harvesting costs. The foremost trait 
that distinguishes popbean from all other types of bean 
is the ability to expand the cotyledons after grains ex- 
plode in response to heating, which is referred to as 
popping expansion, similar to popcorn, although the 
popping mechanism is different. In popcorn, the endo- 
sperm is liquefied and explosive pressure builds up in 
the pericarp [6,7], and both pericarp thickness and endo- 
sperm starch type have been attributed to the popping 
ability. Thus, there are two crucial factors that influence 
the popping of popcorn: whether a kernel can pop or 
not, and if so, to what extent. Popping rate and flake size 
seemingly correspond to these two factors, which to- 
gether determine high popping volume [8]. In contrast, 
popping in nuna beans is the result of pressurized steam 
trapped within and between the mesophyll cells in the 
cotyledons [9,10]. The nutritive value analysis of nuna 
bean revealed that they have a lower mean content of 
protein, phosphorus, iron, and boron than dry bean var- 
ieties, and a higher level of copper and starch, which 
may be related to their unique texture and taste similar 
to roasted peanuts. Lectins and other anti-nutritional 
compounds were higher in raw and boiled nuna samples 
than in toasted nunas, while tannin levels did not change 



from raw to toasted treatments. Overall in-vitro digest- 
ibility was slightly lower for toasted nunas than boiled 
dry bean [11]. Taken together, the nutritional features 
of toasted nuna beans make them a healthy snack, 
although commercial production would require the 
genetic improvement of other agronomic traits, particu- 
larly the day-length sensitivity that has likely restricted 
production and commercialization of nuna beans in 
temperate regions [12,13]. 

In maize, experimental evidence indicate that popping 
characteristics are quantitatively inherited [14,15], con- 
trolled by multiple genes [15-17], and influenced by en- 
vironmental effects [18] and popping methods [19]. Both 
additive and dominant genetic effects play very import- 
ant roles in the inheritance of popping characteristics 
[14], and several putative quantitative trait loci (QTLs) 
have been identified [8,20,21]. However, most of these 
studies on popping ability have focused on the maize 
crop, and there is little information on the ability to con- 
tribute towards popping expansion of the nuna bean. 
The identification of genomic regions associated with 
this popping ability would enable breeders to develop 
improved cultivars using marker assisted selection 
(MAS). To identify these genomic regions, it is import- 
ant not only to establish accurate phenotyping methods, 
but also to develop a saturated molecular marker-based 
genetic linkage map, and then to detect QTLs associated 
with these popping traits. Molecular markers have 
emerged as powerful tools not only for mapping genes/ 
QTLs governing economically important traits in crops 
[22], but also for unlocking the useful genetic diversity 
from unadapted/wild/unrelated germplasm [23]. 

Common bean is a diploid species with a genome size 
estimated at 450 to 650 Mb [24] that is distributed 
among 22 chromosomes (n = 11). The first core genetic 
linkage map of common bean was based on a recombin- 
ant inbred line (RIL) population resulting from the cross 
between representatives of the Mesoamerican (BAT93) 
and the Andean (Jalo EEP558) gene pools, which 
included 194 restriction fragment length polymorphic 
(RFLP) markers [25]. Nowadays the availability of micro- 
satellite markers and their potential for anchoring new 
genetic maps have allowed a new expanded version of 
the core linkage map to be created with several hundred 
of these types of markers, including markers with puta- 
tive gene functions. Thus, the current core common 
bean linkage map covers a genetic distance of 1258.8 cM 
and includes a total of 413 loci placed across 11 linkage 
groups (LGs) with an average distance between neigh- 
boring loci of 3.0 cM [26]. In addition, about fifteen RIL 
mapping populations and more than twenty five link- 
age maps have been developed, most of them created 
from inter-gene pool crosses, which include divergent 
parents showing high genetic polymorphism [27,28]. 
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Nevertheless, few populations have been generated from 
intra-gene pool crosses [29-35], likely due to their low 
polymorphism level [36,37] although this kind of popu- 
lation does not usually show phenotypic abnormalities 
or undesirable segregating individuals [38,39]. Recently, 
a Mesoamerican saturated intra-gene pool map has been 
constructed by combining the genetic information from 
intra- and inter-gene pool segregating populations [35] . 

Mapping QTLs is of great importance to understand 
the genetic architecture underlying complex traits. Pop- 
ping expansion seems to involve several genes, environ- 
mental factors, and gene-gene and gene-environment 
interactions. Therefore, to better understand the genetic 
control of popping expansion, in addition to detecting 
single-locus QTL effects, there is a need to identify the 
interactions among the different loci (epistatic QTLs, E- 
QTLs), as well as their environment interaction effects 
(QTLs x Environment, QE; and E-QTLs x Environment, 
E-QE). Indeed, these interactions have been successfully 
analysed for other complex traits in several crops species 
e.g. rice [40], wheat [41], cotton [42] or maize [43]. 

In order to identify genetic determinants of popping 
ability traits in common bean, we have evaluated the 
popping dimension index (PDI), expansion coefficient 
(EC) and percentage of unpopped seeds (PUS). The first 
two traits reflect structural changes of seeds when heat- 
ing (cotyledon expansion), while the third is of consider- 
able importance for breeding and from a commercial 
point of view. Besides, we have constructed a genetic 
map from an Andean RIL population, which includes 



193 molecular markers. This has permitted the location 
of single-locus and epistatic QTLs involved in PDI, EC 
and PUS, some of which co-localized in four LGs. Mar- 
kers associated to these QTLs could be used as genetic 
tools for MAS programs devoted to popping improve- 
ment of nuna bean cultivars. 

Results 

Phenotypic variation of popping traits and correlations in 
RIL population 

Mean values, standard deviations and ranges of variation 
for popping traits are shown for each environment in 
Table 1. PDI, EC and PUS were clearly different between 
both parents and varied significantly among RILs. In 
addition, these popping traits showed a pattern of con- 
tinuous distribution, which weakly departed from normal- 
ity (data not shown). Strong transgressive segregations 
were observed for these traits, indicating that alleles with 
positive effects are distributed among the parents. The 
broad sense heritability value of EC was rather low (0.33 ± 
0.04), whereas the heritability values of PDI (0.53 ± 0.04) 
and PUS (0.48 ± 0.04) were moderate, suggesting that gen- 
etic gains could be obtained when selecting for these traits. 
The phenotypic correlations between PUS and two pop- 
ping traits showed that PUS was significant and negatively 
correlated with PDI and EC, while EC was significant and 
positively correlated with PDI in all environments assayed 
(Table 2). Highly significant differences were found among 
genotypes and genotype x environment interactions for all 
popping traits (Table 3). 



Table 1 Mean values of the popping traits analysed in the RIL population PMB0225 x PHA1037 



Trait 


Environment 




Parents 








RILs 




PMB0225 


PHA1037 


pa 

r par 


Mean b 


Range 


Pril 


Popping dimension index (PDI) 


LD09 


-0.99 ± 3.28 






2.45 


± 7.69 


-16.21 -45.08 






SD09 


-0.05 ± 1 .44 


22.29 ± 6.18 




3.53 


± 6.17 


-20.83 - 26.37 






LD10 


-0.60 ± 0.30 






2.77 


± 7.30 


-21.95 - 39.91 






SD10 


-1.17 ± 5.37 


27.84 ± 4.81 




1.21 


± 5.53 


-10.44 - 28.96 






LD11 


1.01 ± 1.35 






3.28 


± 7.75 


-26.57 - 31.88 




Expansion coefficient (EC) 


LD09 


2.25 ± 2.63 






11.82 


± 19.22 


0.00 - 95.00 


X* 




SD09 


7.15 ± 10.12 


68.70 ± 22.23 




24.89 


± 37.87 


-53.33 - 383.33 






LD10 


10.19 ± 9.17 






13.37 


± 21.38 


0.00 - 166.67 






SD10 


3.33 ± 4.35 


45.83 ± 9.48 




8.96 


± 16.08 


-28.21 - 105.00 






LD11 


1 7.78 ± 1 .92 






19.65 


± 24.87 


0.00 - 150.00 




Percentage of unpopped seeds (PUS) 


LD09 


95.04 ±4.18 






84.41 


± 25.58 


0.00 - 100.00 






SD09 


95.00 ± 7.07 


24.44 ± 19.44 




52.85 


± 31.78 


0.00 - 100.00 






LD10 


77.78 ± 6.29 






73.12 


± 28.19 


0.00 - 100.00 


-X- 




SD10 


100.00 ±0.00 


14.00 ±0.00 




86.93 


± 27.18 


0.00 - 100.00 






LD11 


97.00 ± 5.24 






77.42 


± 29.99 


0.00 - 100.00 





a * ** Significant differences between the two parents (par.) or among lines (RILs) at P <0.05 and P <0.001, respectively. 
b Mean of 2 and 3 replicates ± standard deviation. 

No data taken for popping traits in the parent PHA1037 under long-day conditions. 
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Table 2 Phenotypic correlations among the popping 
traits evaluated in the RIL population PMB0225 x 
PHA1037 



Trait 


Environment 


PDI 


EC 


EC 


LD09 


0.76** 






SD09 


0.53** 






LD10 


0.75** 






SD10 


0.76** 






LD11 


0.79** 




PUS 


LD09 


-0.80** 


-0.75** 




SD09 


-0.68** 


-0.44** 




LD10 


-0.64** 


-0.53** 




SD10 


-0.86** 


-0.78** 




LD11 


-0.85** 


-0.84** 



*, ** Significant correlation at P <0.05 and P <0.001, respectively. 

PDI, popping dimension index; EC, expansion coefficient; PUS, percentage of 

unpopped seeds. 



Genetic map construction 

We screened PMB0225 and PHA1037 parents for DNA 
polymorphism by using amplified fragment length poly- 
morphic (AFLP), simple sequence repeat (SSR) and sin- 
gle nucleotide polymorphism (SNP) markers. A total of 
92 AFLP primer combinations were assayed, which 
allowed the amplification over 3,700 AFLP fragments, 
279 (7.5%) of which were polymorphic between the par- 
ents. According to their inheritance pattern and reliabil- 
ity features, 18 combinations, producing 94 polymorphic 
loci, were selected for genotyping the RIL population. 
Furthermore, 1035 SSR and 251 SNP markers were 
screened, which rendered polymorphism rates of 10.2% 
and 7.2%, respectively. Thus, 106 SSR and 18 SNP poly- 
morphic loci were also analysed in the RIL population 
for map construction. Some of the polymorphic markers 
(9 AFLPs, 11 SSRs and 5 SNPs) were not linked with 
any other marker in the existing map and so they could 
not be mapped. Finally, the genetic map developed from 
the cross PMB0225 x PHA1037 (Figure 1) was con- 
structed with a total of 193 loci (85 AFLP, 95 SSR, and 
13 SNP markers), of which 101 were dominant and 92 



codominant, resulting in the formation of 12 LGs. These 
LGs were designated according to Pedrosa-Harand et al. 
[44] on the basis of 55 previously mapped common SSR 
markers [26,31-33,45-50], which were used as a guide 
for the assignment of LG number and orientation. The 
map spanned a total genetic distance of 822.1 cM, with 
an average of 68.5 cM per LG, ranging from 16.5 cM 
(LG 6) to 106.4 cM (LG 3). The density of markers ran- 
ged from 1.3 cM (LG 6) to 6.75 cM (LG 7), with an aver- 
age of 4.3 cM per LG. A detailed description of this map 
is provided in Figure 1 and Table 4. 

Single environment QTL analysis of popping traits 

QTL analysis based on MQM mapping using MapQTL 
was undertaken to identify single-locus QTLs in individ- 
ual environments; thus, a total of sixteen QTLs for pop- 
ping traits were detected (Figure 1). Five single-locus 
QTLs were identified for PDI, of which one (PDI3 PP ) 
was identified in four environments, two were detected 
in two environments (PDI5 PP and PDI7 PP ), and the other 
two were detected in only one environment (PDI1 PP and 
PDI6 PP ). The percentage of phenotypic variation 
explained by a single QTL identified for PDI ranged 
from 7.5 (for PDI3 PP in LD09) to 18.7% (for PDI6 PP in 
SD10). For EC, two QTLs were detected in two or more 
environments (EC5 PP and EC7 PP ) and two QTLs were 
identified in only one environment (EC3 PP and EC6 PP ); 
in SD09 no significant QTLs were detected. The per- 
centage of phenotypic variance for EC ranged from 7.3 
(for EC7 PP in LD09) to 15.9% (for EC6 PP in SD10). In 
the case of PUS, seven single-locus QTLs were identi- 
fied, four of which were detected in two environment 
conditions (PUS3 PP , PUS5 PP , PUS6 PP , and PUS7 PP ) and 
three in only one environment (PUS1.2 PP , PUS9 PP , and 
PUS10 PP ). The percentage of phenotypic variance for 
PUS ranged from 6.3 (for PUS1.2 PP and PUS10 PP , both 
detected in LD11) to 16.5% (for PUS6 PP in SD10). For 
each popping trait, QTLs with positive (alleles from 
PHA1037) and negative (alleles from PMB0225) additive 
values were identified. A complete description of the 
MQM mapping analysis is provided in Table 5. 



Table 3 ANOVA for the popping traits measured in the RIL population PMB0225 x PHA1037 


Source of variation 


Degree of freedom 


PDI 


EC 


PUS 


Environment 


4 


367.95** 


8648.06* 


28427.43** 


Block (Environment) 


4 


44.31* 


901.49 


1 703.72** 


Genotype 


186 


188.53** 


1291.77** 


2406.78** 


Genotype x Environment 


602 


41 .80** 


593.09** 


562.39** 


Error 


865 


21.33 


387.05 


392.90 



*, ** Significant differences at P < 0.05 and P < 0.001, respectively. 

PDI, popping dimension index; EC, expansion coefficient; PUS, percentage of unpopped seeds. 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 Genetic linkage map of common bean based on the RIL population PMB0225 x PHA1037. Location of single-locus QTLs and 
E-QTLs controlling popping traits: popping dimension index (PDI), expansion coefficient (EC), and percentage of unpopped seeds (PUS). Common 
SSR markers to previously published maps (see text for references) are indicated in bold. Names of QTLs are listed in Tables 5, 6, and 7. 
Cumulative distances among markers are indicated in cM to the left of the linkage group, names of markers are shown on the right. QTLs are 
depicted as vertical bars to the right of the linkage groups. QTLs detected by both MapQTL and QTLNetwork software packages are indicated in 
black, QTLs identified only by MapQTL are shown in red, and QTLs detected only by QTLNetwork are represented in blue. Epistatic interactions 
between QTLs are represented with numbered stars. 



Multiple environment QTL analysis of popping traits 

In addition to the single-locus QTLs identified for each 
environment, a single-locus QTL analysis for multi- 
environment was also carried out by using QTLNet- 
work. As the result of one-dimensional MCIM analysis, 
a total of fourteen single-locus QTLs were detected 
with significant genetic main effects and/or QTLs x En- 
vironment (QE) effects (Figure 1). Four single-locus 
QTLs were detected for PDI, three of which showed 
only individual additive effects (PDI3 PP , PDI5 PP , and 
PDI6 PP ); and the remaining QTL (PDI7 PP ) had both in- 
dividual additive and QE interaction effects. The per- 
centage of phenotypic variation explained for PDI by a 
single QTL ranged from 2.2 (for PDI6 PP ) to 7% (for 
PDI3 PP ); in the case of QTL PDI7 PP , the phenotypic 
variation explained by QE interaction effects ranged 
from 0.6 (in SD09) to 1.6% (in SD10). For EC, four sin- 
gle locus QTLs were detected, two of which displayed 
only individual additive effects (EC3 PP and EC5 PP ) while 
the other two showed both individual additive and QE 
interaction effects (EC7 PP and EC9 PP ). The percentage 
of phenotypic variance for EC ranged from 2.8 (for 
EC7 PP ) to 4.6% (for EC5 PP ); the phenotypic variation 
explained by QE interaction effects ranged from 0.8 (for 



EC9 in SD10) to 1.4% (for EC7 in SD10). Six QTLs 
were identified for PUS, two of which had only individ- 
ual additive effects (PUS3 PP and PUS5 PP ) while the 
other four displayed both individual additive and QE 
interaction effects (PUS1.1 PP , PUS1.3 PP , PUS6 PP , and 
EC7 PP ). The percentage of phenotypic variance for PUS 
ranged from 0.6 (for PUS1.3 PP ) to 7.4% (for PUS3 PP ); 
the phenotypic variation explained by QE interaction 
effects ranged from 0.8 (for PUS1.1 PP in SD10) to 1.8% 
(for PUS1.3 PP in SD10). For the three popping traits, 
single-locus QTLs with positive and negative additive 
values were detected, indicating that alleles from both 
parents, PHA1037 and PMB0225, have a positive agro- 
nomical effect on the traits. A complete report of the 
one-dimensional genome scan analysis using QTLNet- 
work is given in Table 6. 

Epistatic and environmental interactions 

Given that popping expansion seem to be a complex poly- 
genic trait, genetic interactions may have significant effects 
on the phenotypic values. Therefore, two-dimensional gen- 
ome scan was undertaken for multi-environment analysis 
using QTLNetwork to identify epistatic and environment 
interactions among QTLs. A total of ten epistatic QTLs 



Table 4 Distribution of molecular markers on the linkage map constructed from the RIL population PMB0225 x 
PHA1037 



Linkage groups 


Map length (cM) 


No. of markers 


Marker density (cM/marker) 


AFLP 


Marker types 
SSR 


SNP 


1a 


95.6 


15 


6.4 


9 


6 




1b 


60.7 


19 


3.2 


7 


10 


2 


2 


64.5 


20 


3.2 


7 


11 


2 


3 


106.4 


23 


4.6 


7 


16 




4 


64.9 


14 


4.6 


3 


5 


6 


5 


46.1 


8 


5.8 


3 


5 




6 


16.5 


13 


1.3 


11 


2 




7 


81.1 


12 


6.8 


6 


6 




8 


77.6 


20 


3.9 


8 


12 




9 


90.3 


14 


6.5 


5 


9 




10 


54.5 


15 


3.6 


8 


6 


1 


11 


63.9 


20 


3.2 


11 


7 


2 


Total 


822.1 


193 


4.3 


85 


95 


13 
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Table 5 Single-locus QTLs detected for popping traits using multiple-QTL model mapping for individual environment 
analysis 



Trait Environment 


QTL 

(position) a 


Linkage 
group 


Closest marker 
(position) b 


LOD score c 


LOD threshold d 


R 2e 


A f 


Popping dimension index (PDI) 
















LD09 


PDI3 PP (68.1-74.5) 


3 


IAC24 (71.0) 


2.7 


2.6 


7.5 


1.92 




PDI5 PP (32.8-38.4) 


5 


BMc321 (37.8) 


4.3 




13.4 


2.61 


SD09 


PDI1 PP (55.6-59.3) 


1b 


E32M60-147 (58.2) 


3.0 


2.8 


7.7 


1.57 




PDI3 PP (62.8-66.3) 


3 


BMd1 (63.7) 


4.4 




12.6 


1.97 


LD10 


PDI3 PP (58.7-65.8) 


3 


PVEST042 (60.0) 


3.1 


2.6 


8.2 


1.93 




PDI7 PP (33.4-37.8) 


7 


BMc294 (36.6) 


4.7 




12.9 


-2.44 


SD10 


PDI3 PP (62.1-65.9) 


3 


BMd1 (63.7) 


4.7 


2.7 


9.2 


1.46 




PDI6 PP (1.8-5.5) 


6 


E40M60-91 (5.0) 


8.9 




18.7 


2.15 


LD11 


PDI5 PP (35.9-40.4) 


5 


BM138 (39.8) 


3.0 


2.8 


7.7 


2.07 




PDI7 PP (22.9-37.6) 


7 


BMc294 (36.6) 


5.1 




13.3 


-2.72 


Expansion coefficient (EC) 
















LD09 


EC5 PP (41.3-46.1) 


5 


BM175 (46.1) 


3.6 


2.6 


9.8 


6.18 




EC7 PP (30.1-37.7) 


7 


BMc294 (36.6) 


2.7 




7.3 


-5.25 


SD09 


Not significant 








2.1 






LD10 


EC7 PP (30.5-37.9) 


7 


BMc294 (36.6) 


3.5 


2.4 


10.7 


-6.62 


SD10 


EC3 PP (67.2-76.8) 


3 


IAC24 (71.0) 


5.9 


2.6 


12.0 


5.69 




EC6 PP (1.2-5.4) 


6 


E40M60-91 (5.0) 


7.6 




15.9 


6.75 


LD11 


EC5 PP (39.1-43.3) 


5 


BM138 (39.8) 


4.2 


2.8 


9.8 


6.82 




EC7 PP (32.7-41.6) 


7 


BMc294 (36.6) 


5.9 




14.2 


-8.18 


Percentage of unpopped seeds (PUS) 
















LD09 


PUS5 PP (38.3-46.1) 


5 


BM175 (46.1) 


3.9 


2.9 


12.3 


-9.14 


SD09 


PUS3 PP (66.9-74.6) 


3 


PvM152a (67.6) 


3.3 


2.9 


8.9 


-9.11 




PUS9 PP (42.3-46.1) 


9 


E31M51-59 (27.5) 


3.4 




10.1 


9.78 


LD10 


PUS6 PP (3.1-6.4) 


6 


E40M60-91 (5.0) 


3.2 


2.9 


9.7 


-8.63 




PUS7 PP (22.7-38.4) 


7 


BM185 (24.2) 


3.0 




8.9 


7.81 


SD10 


PUS3 PP (67.1-70.6) 


3 


PvM152a (67.6) 


4.5 


2.8 


8.8 


-8.26 




PUS6 PP (1.6-6.2) 


6 


BMc238 (2.3) 


7.7 




16.5 


-11.36 


LD11 


PUS1.2 pp (35.9-39.1) 


1b 


E43M38-138 (36.4) 


3.3 


2.8 


6.3 


-7.23 




PUS5 PP (38.1-40.2) 


5 


BM138 (39.8) 


4.6 




8.7 


-8.33 




PUS7 PP (26.9-39.3) 


7 


BMc294 (36.6) 


5.2 




9.9 


9.29 




PUS10 PP (6.2-7.1) 


10 


BMb414 (7.0) 


3.4 




6.3 


-6.79 



Those QTLs identified in 2 or more environments are marked in bold. 

a QTL name (according to Miklas and Porch 2010) and its estimated map position (in Kosambi cM). 
b Nearest marker to peak of the detected QTL and its map position (in Kosambi cM). 
c LOD score detected for the nearest marker. 

d LOD thresholds score determined by permutation test for each trait in each environment (P = 0.05). 
e Percentage of the total phenotypic variation explained by the QTL. 

f Estimated additive effect. Positive values indicate allele arising from PHA1037 and negative values indicate allele arising from PMB0225. 



(E-QTLs) involved in six epistatic interactions (Table 7 
and Figure 1) were detected for the three evaluated 
traits. The percentage of phenotypic variance explained 
by the interaction of these E-QTLs was low, ranging 
from 0.2 (for EC) to 2.7% (for PDI). Interestingly, the 
ten E-QTLs identified were previously detected as 
single-locus QTLs, which indicated that these QTLs not 



only participated in epistatic interactions, but they also 
had an individual effect. Regarding PDI, only one epi- 
static interaction was detected, between QTLs E-PDI3 PP 
and E-PDI7 PP , explaining 2.7% of phenotypic variation. 
For EC, two epistatic interactions were identified in- 
volving four E-QTLs; the percentage of phenotypic vari- 
ance explained by the interactions of E-EC3 PP -E-EC9 PP 
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Table 6 Single-locus QTLs and QTLs x Environment (QE) effects detected for the popping traits using multi- 
environment analysis 



QTL 



Marker 
interval 



LG 

(position) 3 



/7 2 (a) c 



QE AE 



h 2 (aef 



Popping dimension index (PDI) 

PDI3 PP 

PDI5 PP 

PDI6 PP 

PDI7 PP 



Expansion coefficient (EC) 

EC3 PP 

EC5 PP 

EC7 PP 



EC9 PP 

Percentage of unpopped seeds (PUS) 

PUS1.1 PP 

PUS1.3 PP 

PUS3 PP 
PUS5 PP 
PUS6 PP 

PUS7 PP 



IAC24-BM287 
E32M60-263-BMc321 
BMc238-E40M60-91 
BM185-BMc294 



BMc259-IAC24 

E32M60-100-BM175 

BM185-BMc294 

PV-at007-BMc184 

BMc324-BM200 
E32M60-147-SNP-5503 

BMc259-IAC24 

BM138-E32M60-100 

E40M60-91-E45M50-50 

BMc294-BMc248 



3 (71.0-84.9) 

5 (35.8-37.8) 

6 (2.3-5.0) 

7 (24.2-36.6) 



3 (67.8-71.0) 
5 (40.6-46.1) 
7 (24.2-36.6) 

9 (60.9-70.2) 

1a (66.5-95.6) 
1 b (58.2-60.6) 

3 (67.8-71.0) 

5 (39.8-40.6) 

6 (5.0-5.8) 

7 (36.6-38.1) 



3.74 
1.31" 
1 .46** 
-1.75" 



3.67 
4.08** 
-6.14*' 

3.39** 

1 0.29*' 
-7.01" 

-11.51" 
-8.04** 



4.99 



7.0 
6.9 
2.2 
6.1 



4.3 
4.6 
2.8 

3.4 

1.6 
0.6 

7.4 
6.5 
4.8 

3.1 



ns 

1.01* AE2 
1.41" AE4 
-1.44** AE5 



ns 
ns 

5.77** AE4 
-4.58* AE5 

-2.75* AE4 



-10.92 AE4 

-8.29** AE2 
7.39** AE4 

ns 

ns 

-3.85** AE4 
3.58* AE5 

-3.17* AE4 
3.51* AE5 



0.6 
1.6 
1.4 



1.4 
0.9 

0.8 



0.8 



1.8 



1.6 
1.1 

1.2 
1.1 



a Linkage group and the estimated confidence interval of QTL position in brackets (in Kosambi cM). 

b Estimated additive effect. Positive values indicate that alleles from PHA1037 have a positive effect on the traits, and negative values indicate that positive effect 

on the traits is due to the presence of the alleles from PMB0225. 

c Percentage of the phenotypic variation explained by additive effects. 

d Predicted additive by environment interaction effect. AE1, AE2, AE3, AE4, and AE5 additive by environment interaction effect associated with environments 
LD09, SD09, LD10, SD10, and LD1 1, respectively. The meaning of sign values is described in the second footnote ( b ). 
f Percentage of the phenotypic variation explained by additive x environment interaction effect. 

*P <0.05, **P <0.01, ***P < 0.001. Only significant effects are listed, ns = No significant effects on the five environmental conditions evaluated. 



and E-EC5 PP -E-EC7 PP was 0.2% and 2.6%, respectively. E- 
QTL x Environment (E-QE) interaction effects were not 
significant for either PDI or EC popping traits. In the case 
of PUS, three epistatic interactions and four E-QTLs were 
detected, and two of the interactions showed only genetic 
effects and explained 0.6% (E-PUS3 PP -E-PUS7 PP ) and 
1.8% (E-PUS5 PP -E-PUS6 PP ) of phenotypic variation, re- 
spectively. The remaining epistatic interaction (E-PUS3 PP - 
E-PUS6 PP ) had both genetic and E-QE interaction effects, 
and the percentage of phenotypic variation obtained was 
1.1% and 0.9% for genetic and E-QE interaction effects, re- 
spectively. A detailed description of the digenic epistatic 
interaction analysis is shown in Table 7. 

Discussion 

The main goal of the current study was to unravel the 
genetic architecture of popping ability in nuna bean. 
Thus, popping traits related to changes in the physical 
structure of seeds have been analysed on the basis of 



their similarity to popcorn, whose cotyledons also expand 
when dry grains are heated. Genetic analysis performed 
indicated that popping ability traits show a polygenic in- 
heritance, making this the first work to report the genetic 
control of these traits in common bean. Transgressive seg- 
regation was observed for popping traits, suggesting that 
combinations of alleles from both parents have effects in 
the same direction; in fact, not only PHA1037 but also 
PMB0225 bear alleles with a positive effect on popping 
ability, a finding backed up by QTL analyses. Since trans- 
gressive segregation relies on additive genetic variation, 
the extreme phenotypes can be maintained and fixed 
through artificial selection, providing the potential for im- 
provement of popping ability. Furthermore, the analysis of 
variance showed that although the genotype x environ- 
ment interaction affects popping ability, this effect is fairly 
uniform across all genotypes, and it does not seriously 
compromise genotypic main effects, making progress 
from selection feasible. 



Table 7 Epistatic QTLs (E-QTLs) and E-QTL x Environment (E-QE) effects detected for popping traits using multi-environment analysis 

E-QTLi a Marker interval LG (position) 13 E-QTLj a Marker interval LG (position) AA C /? 2 (aa) d E-QE AAE e /? 2 (aae) f 

Popping dimension index (PDI) 



E-PDI3 PP 


IAC24 -BM187 


3 (71.0-84.9) 


E-PDI7 PP 


BM185 -BMc294 


7 (24.2-36.6) 


-2.98*** 


2.7 


ns 


Expansion coefficient (EC) 


















E-EC3 PP 


BMc259 -IAC24 


3 (67.8-71.0) 


E-EC9 PP 


PV-at007 -BMc184 


9 (60.9-70.2) 


1.62* 


0.2 


ns 


E-EC5 PP 


E32M60-100 -BM175 


5 (40.6-46.1) 


E-EC7 PP 


BM185 -BMc294 


7 (71.0-84.9) 


-7.59*** 


2.6 


ns 


Percentage of unpopped seeds (PUS) 


















E-PUS3 PP 


BMc259 -IAC24 


3 (67.8-71.0) 


E-PUS6 PP 


E40M60-91 -E45M50-50 


6 (5.0-5.8) 


-3.87*** 


1.1 


-4.61* AAE4 0.9 


E-PUS3 PP 


BMc259 -IAC24 


3 (67.8-71.0) 


E-PUS7 PP 


BMc294 -BMc248 


7 (36.6-38.1) 


5.47*** 


0.6 


ns 


E-PUS5 PP 


BM138 -E32M60-100 


5 (39.8-40.6) 


E-PUS6 PP 


E40M60-91 -E45M50-50 


6 (5.0-5.8) 


-3.62*** 


1.8 


ns 



a E-QTLi and E-QTLj are the two QTLs involved in epistatic interaction. 

b Linkage group and the estimated confidence interval of QTL position in brackets (in Kosambi cM). 

c Estimated additive by additive epistatic effect. Positive values indicate that alleles from PHA1037 have a positive effect on the traits, and negative values indicate that positive effect on the traits is due to the 
presence of the alleles from PMB0225. 

d Percentage of the phenotypic variation explained by additive by additive epistatic effects. 

e Predicted additive by additive epistatic effect by environment interaction effect. AAE1, AAE2, AAE3, AAE4, and AAE5: epistasis associated with environments LD09, SD09, LD10, SD10, and LD1 1, respectively. The 
meaning of sign values is described in the third footnote ( c ). 

f Percentage of the phenotypic variation explained by additive by additive epistatic effect by environment interaction effect. 

*P <0.05, **P <0.01, ***P < 0.001. Only significant effects are listed, ns = No significant effects on the five environmental conditions evaluated. 
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A comprehensive QTL analysis was performed to de- 
tect single-locus QTLs, epistatic QTLs and their envir- 
onment interactions on a newly created genetic linkage 
map. This map was constructed for an Andean intra- 
gene pool cross involving PMB0225 (dry bean) and 
PHA1037 (popbean) parents. Despite the morphological 
diversity observed in the Andean intra-gene pool, the 
low genetic polymorphism existing in this common bean 
germplasm hinders the development of genetic linkage 
maps [31,33,51,52]. Our results confirmed a low poly- 
morphism between the Andean parents, thus the overall 
polymorphism rate detected was 8.3%. Likewise, Blair 
et al. [33] screened a total of 700 SSR markers on the 
Andean parents G21242 and G21078, but only 74 map- 
pable markers were found in that survey resulting in a 
polymorphism rate of 10.6%, comparable to the level of 
polymorphism here reported (i.e. 10.2% for SSR markers). 
However, Cichy et al. [31] found a moderate SSR poly- 
morphism of 30% between the Andean parents used to 
generate a G 19833 x AND696 RIL mapping population. 

Interestingly, the genetic map described in this work 
shares 55 SSR markers with previously published com- 
mon bean maps [26,31-33,45-50]. In fact, linkage asso- 
ciations have been found in terms of SSR marker 
mapping in the present map and previous maps, while 
the collinear order of the commonly mapped SSR loci 
has been generally observed although some inversions 
affecting SSR markers located close to one another 
have been detected. Since these loci were distributed 
throughout all LGs, they would permit the alignment 
of homologous LGs between maps and facilitate mar- 
ker transfer across populations as well as between 
related species. Hence, these shared markers could be 
used as anchor points for map merging and syntenic 
analysis such as Galeano et al. [35] have recently 
reported for the consensus Mesoamerican intra-gene 
pool map. 

The genetic linkage map developed herein includes 
193 loci (85 AFLP, 95 SSR, and 13 SNP markers) across 
12 LGs that cover a genetic distance of 822.1 cM, with 
an average of 4.3 cM per marker. Prior to this work, two 
Andean maps have been described for QTL analysis; the 
map depicted by Cichy et al. [31] included 167 loci that 
spanned a total map length of 1105 cM with an average 
marker density of 6.6 cM per locus. On the other hand, 
the map constructed by Blair et al. [33] contained 118 
loci with a total map length of 726.0 cM and a mean 
marker density of 6.2 cM per locus. Therefore, com- 
pared to previous Andean maps [31,33], the genetic map 
here developed shows a suitable marker density and gen- 
ome coverage, which has permitted the first identifica- 
tion of popping ability QTLs. 

Three closely related popping traits such as PDI, EC, 
and PUS have been analysed and the reliability of the 



QTLs associated to these traits has been enhanced by 
using several software programs, which decreased 
the risk of detecting false positive and negative QTLs 
[53-55]. Therefore, single and multi- environment QTL 
analyses were performed to dissect the genetic archi- 
tecture of popping ability in nuna bean. In summary, a 
total of nineteen single-locus QTLs were identified by 
MapQTL and QTLNetwork. Eleven of the fourteen QTLs 
identified by QTLNetwork were also detected by 
MapQTL; thus, the two independent approaches con- 
verged on the identification of common single-locus QTLs 
for PDI (PDI3 PP , PDI5 PP , PDI6 PP , and PDI7 PP ), EC (EC3 PP , 
EC5 PP , and EC7 PP ), and PUS (PUS3 PP , PUS5 PP , PUS6 PP , and 
PUS7 PP ). The percentage of phenotypic variation explained 
by the single-locus QTLs identified by MapQTL for pop- 
ping traits was comparatively higher than that of the 
QTLNetwork. The results of multi-environment analyses 
showed that genetic main effects were sometimes subject 
to environmental modification; this could explain why we 
obtained a lower phenotypic variance using QTLNetwork 
software. Even so, for multi-environment analyses, the per- 
centage of phenotypic variance attributable to genetic 
effects was as expected for a complex trait, which is gov- 
erned by several small effect QTLs/genes located in differ- 
ent genomic regions, and where the environment 
interactions play an important role. Therefore, the four 
common single-locus QTLs detected for PDI and PUS to- 
gether explained 22.2 and 21.8% of the phenotypic vari- 
ance, respectively. Regarding EC, the three common QTLs 
explained 11.7% of the phenotypic variance in the RIL 
population. In addition, it was interesting to find that the 
common QTLs (detected with both software programs) 
not only were consistent over environments, but they also 
co-localized with QTLs for the analysed popping traits. 

Overall, significant positive correlations between PDI 
and EC and negative correlations among PUS, and PDI 
and EC, together with the detection of co-localized 
QTLs for PDI, EC, and PUS on LGs 3, 5, 6, and 7, sug- 
gested that QTLs for popping ability are not evenly dis- 
persed throughout the genome but rather are clustered 
in several genomic regions. The QTLs sign values of 
additive effects corresponded to the significant genotypic 
correlations observed among the three popping traits. 
Thus, the co-localized QTLs located on LGs 3, 5 and 6 
showed positive (alleles from PHA1037) and negative 
(alleles from PMB0225) values of additive effects for PDI 
and EC, and PUS, respectively. Meanwhile, the opposite 
sign values of additive effects were found for the co- 
localized QTLs located on LG 7, which indicated that 
PMB0225 also contributes positively to popping ability. 
To date, research into popping ability has focused on 
popcorn. Thus, Babu et al. [21] detected four QTLs for 
popping expansion volume, five for flake volume, and 
five for percentage of unpopped kernels, and revealed 
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QTLs in overlapping or mostly adjoining regions in the 
same chromosomes affecting two or three popping 
traits. Likewise, Li et al [8] evaluated three important 
traits for popcorn (i.e. popping volume, flake size, and 
popping rate), and six chromosome regions were found 
to control two or three popping traits simultaneously. 
Hence, as in nuna bean, the detection of co-localized 
QTLs for popping traits suggested that either pleiotropic 
QTLs controlled several popping traits, or tightly linked 
QTLs for different traits are present together in the 
same genomic regions. The issue of pleiotropy versus 
tight linkage of QTLs may be resolved in the future 
through fine mapping of the target genomic regions. 

Epistatic effects are often involved in complex traits, 
but they are difficult to confirm because of their usually 
small effects and environmental interactions. Genetic 
models for QTL mapping assuming no epistasis could 
lead to biased estimation of QTL parameters, and subse- 
quently result in considerable loss of response in MAS. 
In popcorn, Li et al. [8] carried out a preliminary epi- 
static analysis and detected thirteen pairs of digenic 
interactions for popping ability. In the present work, sev- 
eral epistatic interactions were found involving all of the 
evaluated popping traits. A total of ten E-QTLs, involved 
in six epistatic interactions, were detected, and only one 
epistatic interaction for PUS showed significant E-QE 
interaction effect. Although the phenotypic variation 
explained by each epistatic interaction was found to be 
small, it is interesting to note that the genomic regions 
located on LGs 3, 5, 6 and 7 not only harbor QTLs that 
have individual genetic effects, but are also involved in 
epistatic interactions. Therefore, QTL analysis revealed 
that popping ability of nuna bean is controlled by several 
QTLs, which have only individual additive effects, or 
may also be involved in epistatic or environmental inter- 
actions, indicating that popping is inherited as a poly- 
genic trait, and that epistasis could play a key role. 

Nowadays, popping of common bean is considered an 
interesting agronomic trait, since it allows greater diver- 
sification of this crop as well as the commercialization of 



nuna bean as a new snack product. In popcorn, selection 
for increased expansion coefficient has been successfully 
achieved given its high heritability value [20,21]. Simi- 
larly, Vorwald and Nienhuis [56] estimated that the nar- 
row sense heritability values of fully expanded seeds 
after popping (popping percentage) and expansion coef- 
ficient in nuna bean were relatively high, 0.87 ± 0.07 and 
0.74 ± 0.09, respectively. The broad sense heritability 
values calculated in the present work were moderate, 
suggesting that genetic gain could be obtained for pop- 
ping ability in this legume species. The introgression of 
popping and the development of new day-length insensi- 
tive popbean cultivars would require genetic tools which 
facilitate efficient genotyping selection. Conventional 
phenotype selection methods for popping traits are la- 
borious and time-consuming; consequently, MAS pro- 
vides an efficient and cost-effective alternative that 
accelerates the selection of interesting genotypes. How- 
ever, MAS approaches have been difficult to apply in the 
case of complex traits such as popping ability, because 
individual QTLs have small genetic effects which in 
many cases are also environmentally modulated. Conse- 
quently, the identification of potential candidate QTLs 
for MAS is crucial. Based on the results obtained in our 
study, the co-localized QTLs located on LGs 3, 5, 6, and 
7 are good candidates for MAS, since they showed sta- 
bility across significantly correlated traits, while also 
sharing QTLs for more than one trait, and they could be 
manipulated simultaneously in breeding programs. 
Breeding of nuna cultivars would require adapting them 
to temperate regions, and for this purpose it is import- 
ant to improve their insensitivity to photoperiod. The 
Ppd gene regulates photoperiod sensitivity and it is 
located on LG 1 [57], while popping QTLs are located 
on LGs 3, 5, 6, and 7. Therefore, the use of QTL marker 
assisted selection would facilitate the introgression of 
photoperiod insensitivity without loss of popping ability. 
QTL pyramiding approach would also permit the com- 
bination of QTL alleles with positive effects for popping 
ability on a day-length insensitive genotype through 




Figure 2 Unpopped and popped seeds of the parental genotypes. Both PMB0225 (A) and PHA1037 (B) belong to the Andean gene pool of 
common bean. Scale bar 1 cm. 
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molecular breeding, thus overcoming the main drawback 
for the production and commercialization of nuna beans 
in temperate regions. In this research, some RILs 
showed popping expansion ability and flowered inde- 
pendently of photoperiod conditions. These lines consti- 
tute an interesting breeding goal, and they will hopefully 
allow researchers to isolate the genes and to understand 
the molecular and physiological mechanisms underlying 
agronomic traits which are relevant for the genetic im- 
provement of nuna beans. 

Conclusions 

We have developed a novel Andean genetic linkage map, 
which has permitted the first identification of popping 
ability QTLs in common bean. Our results revealed 
that popping ability of nuna bean is controlled by sev- 
eral QTLs, which have only individual additive effects 
or may also be involved in epistatic or environmental 
interactions, indicating that popping is inherited as a 
polygenic trait, and that epistasis could play a key 
role in its genetic control. Individual and multi- 
environment QTL analyses detected a total of nineteen 
single-locus QTLs, most notably those co-localized for 
the three popping ability traits placed on LGs 3, 5, 6, and 
7. These QTLs showed an individual effect and also 
participated in epistatic interactions. Consequently, the 
co-localized QTLs for popping expansion response are 
useful tools for MAS breeding programs intended to im- 
prove production and adaptation of nuna bean cultivars. 
The results here reported can contribute towards the 
diversification of the nuna bean crop, which is becom- 
ing increasingly relevant as a new food product in the 
agro-food industry due to its nutritional and healthy 
properties. 

Methods 

Population development 

Parents included two genotypes from the Andean gene 
pool, PMB0225 and PHA1037, belonging to Nueva 
Granada and Peru races, respectively. PMB0225 is a 
Spanish improved cultivar resistant to the bean common 
mosaic virus, which shows indeterminate erect growth 
habit type II, white flowers and large seeds. PHA1037 is 
a photoperiod-sensitive nuna popbean germplasm acces- 
sion from Bolivia that has purple flowers and large red 
seeds, and possesses an indeterminate climbing growth 
habit type IV (Figure 2). Both parents were chosen to fa- 
cilitate the introgression of the popping trait in common 
bean cultivar adapted to temperate areas. Individual 
plants of a F 2 segregating population generated from the 
cross of PMB0225 x PHA1037 were selfed to develop 
185 F 7 recombinant inbred lines (RILs) by single-seed 
descent. 



Experimental design 

The 185 RILs and both parents were evaluated for pop- 
ping seed traits using 15 plants per genotype in five 
greenhouse environments over three consecutive years 
(2009, 2010 and 2011). Plants were grown under long- 
day (LD) and short-day (SD) natural photoperiod condi- 
tions with average day and night temperatures of 25°C 
and 20°C, respectively. Sowing dates of LD experiments 
were February 20, 2009 (LD09 code), March 15, 2010 
(LD10 code), and April 27, 2011 (LD11 code), while 
sowing dates of SD experiments were August 15, 2009 
(SD09 code) and September 21, 2010 (SD10 code). For 
all environments the experiments were conducted in a 
randomized complete block design with two or three 
replicates of single row plots (3.0 x 0.8 m). Each plot 
was sown with two seeds per hill and adjusted to a crop 
density of about 30,000 plants/ha. Pods were harvested 
when they were completely dried. Seeds were removed 
and cleaned using a mechanical thrasher followed by 
hand cleaning and winnowing. Then, seeds were stored 
at 5°C for about one month before initiation of the 
present research. 

Trait measurements 

Three popping component traits have been measured: 
popping dimension index (PDI), popping expansion co- 
efficient (EC), and percentage of unpopped seeds 
(PUS). Samples of 50 seeds from each treatment, repli- 
cate and environment, were popped in a Palson 
Denver popcorn maker (1200 W, 230 V, 50 Hz) for 
150 s. Seed was considered fully popped when cotyle- 
dons had expanded sufficiently to shed the seed coat, 
and unpopped or partially popped when the seed coat 
failed to crack or no expansion of the cotyledons was 
observed. Seed dimensions were scored for each indi- 
vidual seed: length (mm) was measured as the longest 
distance across the seed parallel to the hilum, height 
(mm) as the longest distance perpendicular to length, 
and width (mm) as the longest distance across the 
hilum seed. Seed length, width and thickness before and 
after popping were determined from a random 10-seed 
sample, and PDI was recorded as [(^popped seed 
dimensions-^unpopped seed dimensions) /^unpopped 
seed dimensions) x 100]. Each 50-seed sample was 
placed in a graduated cylinder and distilled water was 
added to a total volume of water and seeds of 100 mL. 
The total volume of water added was subtracted from 
the total volume to give the unpopped seed volume 
(unPV). The seeds were drained and patted dry with 
paper toweling and immediately popped to minimize 
absorption of water by the seeds. The volume of seeds 
after popping (PV) was obtained using a procedure 
similar to that used for unpopped seeds. The EC was 
defined as [(PV-unPV)/unPV] x 100. The PUS was 
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calculated as the proportion of 50 seeds that were 
unpopped. 

Statistical data analysis 

Analysis of variance with the Generalized Linear Model 
(GLM) procedure was applied to analyse phenotypic 
data using the SAS Software [58]. Single degree-of- 
freedom orthogonal contrasts between parents were cal- 
culated to show significant differences. For each trait 
and environment, mean value, standard deviation and 
range of variation were calculated. The phenotypic cor- 
relation coefficients between PUS and EC, and PDI, were 
estimated by using PROC CORR [58]. Broad sense herit- 
ability (H 2 ) was estimated as H 2 = o 2 G/(o 2 e/re + o^GE/e + 
o^G), where o^G is the estimate of genotypic or RIL vari- 
ance, o^e is the estimate of error variance, o^GE is the es- 
timate of genotype x environment interaction variance, r 
is the number of replicates per environment, and e is the 
number of environments. The genetic components of vari- 
ance were estimated with the MIXED procedure of SAS 
software [58]. 

DNA extraction and molecular marker analyses 

Total genomic DNA was isolated from young leaves as 
described by Doyle and Doyle [59]. DNA was kept in 
sterile water, visualized after electrophoresis in 1% agar- 
ose gels in IX SB buffer (10 mM sodium boric acid), 
and quantified by comparison with DNA standards 
(Lambda phage DNA digest with Hindlll; Invitrogen Life 
Technologies). DNA was diluted in sterile water to a 
stock concentration of 5-10 ng/[iL and stored at -20°C 
for use in PCR analysis. 

Analysis of AFLP markers was carried out according 
to the procedure described by Vos et al. [60], with some 
modifications [61]. A total of 500 ng of genomic DNA 
was digested with 5 U of Msel and EcoRI enzymes for 
2 h at 37°C in a final volume of 40 [iL. The DNA frag- 
ments were ligated to appropriate adapters via addition 
of 1 U of T4-DNA Ligase (Roche) and incubated for 1 h 
at 37°C. The pre-amplification reactions were performed 
in a volume of 20 \iL using A as the selective nucleotide 
for the Eco primer (Eco + A) and four different Mse pri- 
mers (Mse + A, Mse + C, Mse + G, and Mse + T). The 
PCR cycling parameters were 20 cycles at 94°C for 30 s, 
56°C for 60 s, and 72°C for 60 s. Subsequent PCR ampli- 
fications were performed with primers that included 
three selective bases in their sequences. To detect AFLP 
fragments, Eco primers were labelled using a fluorescent 
dye (FAM, NED, PET or VIC), and the following PCR 
cycling parameters were used for selective amplifica- 
tions: an initial cycle at 94°C for 30 s, 65°C for 30 s, and 
72°C for 60 s. During the next 12 cycles, the annealing 
temperature was lowered by 0.7°C per cycle. The 
temperature conditions for the next 23 cycles were 94°C 



for 30 s, 56°C for 30 s, and 72°C for 60 s. The PCR pro- 
ducts of selective amplifications were separated by capil- 
lary electrophoresis using a DNA sequencer (ABI 
PRISM® 3130 XL Genetic Analyzer, Applied Biosys- 
tems, USA). An internal size marker, GeneScan 500 LIZ 
(35-500 bp; Applied Biosystems) was added, allowing 
the co-loading of different labelled reactions. Data 
regarding selectively amplified DNA fragments were 
analysed with GeneMapper Software 3.7 (Applied Bio- 
systems). Each AFLP marker name included a Keygene 
primer code (http://www.wheat.pw.usda.gov/ggpages/ 
keygeneAFLPs.html) followed by the fragment size in 
base pairs. 

Different sets of SSR and SNP markers previously 
reported (references are shown below) were tested for 
polymorphism in the parental genotypes, and poly- 
morphic loci were used for the construction of the gen- 
etic linkage map. SSR markers were named according to 
the respective authors (IAC-, [49,62,63]; ATA-, [45]; 
BM-, GATS-, [64,65]; BMb-, [47]; BMc-, [46,52]; BMd-, 
[66]; PV-, [67]; PVBR-, [68,69]; PVEST-, [50]; PvM-, FJ-, 
[26,70]). PCR reactions were carried out following the 
protocols described in the publications mentioned 
above, although PCR conditions were changed for some 
SSR markers. Data analysis of the SSR markers was per- 
formed by using either gel electrophoresis or capillary 
electrophoresis in an ABI PRISM® 3130 XL Genetic 
Analyzer (Applied Biosystems, USA). SNP markers 
tested in our mapping population were designated as 
BSNP- [71], Leg- [72], and SNP- [73]. High resolution 
melting technology (HRM) was employed to analyse 
the SNP markers using a LightScanner instrument 
(Idaho Technology), following the protocols described 
by Montgomery et al. [74]. 

Linkage map construction and QTL analyses 

JoinMap® 4.0 software [75] was used to generate the 
linkage maps. Marker data were assigned to LG using a 
minimum logarithm of odds ratio (LOD) score of 6.0, 
and a recombination frequency value of 0.3. The 
Kosambi map function [76] was used to calculate the 
genetic distance between markers. The LGs were desig- 
nated according to Pedrosa-Harand et al. [44]. 

Candidate QTL regions for popping traits were identi- 
fied by using two different mapping software packages, 
MapQTL® 5.0 [77] and QTLNetwork 2.0 [78]. Interval 
mapping and multiple QTL model (MQM) mapping 
were used to detect single-locus QTLs for each environ- 
ment separately by MapQTL. Thus, once potential QTLs 
were detected by interval mapping analysis, markers 
with higher LOD scores were selected as cofactors and 
tested using the automatic cofactor selection procedure 
(default P value cut off for elimination of a cofactor set 
of 0.02). Using the set of selected cofactors, MQM 
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mapping analyses were carried out. A permutation test 
(1000 cycles) was used to determine the LOD threshold 
score at which the QTL was deemed to be present in a 
particular genomic region with a confidence interval of 
95%. In addition, QTLNetwork software was used to 
identify single locus QTLs, epistatic QTLs (E-QTL) and 
their environment interaction effects (QTLs x Environ- 
ment, QE; and E-QTLs x Environment, E-QE) across 
environments. The mixed-model based composite 
interval mapping method (MCIM) was carried out for 
one-dimensional genome scan to detect putative QTLs 
and their environment interactions, and for two- 
dimensional genome scan to identify epistatic effects. An 
experimental-wise significant level of 0.05 was desig- 
nated for candidate interval selection, putative QTL de- 
tection, and QTL effect. Both testing and filtration 
window size were set at 10 cM, with a walk speed of 
1 cM. The critical F-value to declare putative QTLs was 
determined by 1000 permutation test. The effects of 
QTLs and environment interactions were estimated by 
Markov Chain Monte Carlo method [79]. QTLs with 
only genetic effects indicated that these were expressed 
in the same way across environments. In addition, QTLs 
with environment interaction effects suggested that their 
expressions were environmentally dependent. The 
detected QTLs were designated as recommended by 
Miklas and Porch [80]. The genetic map and the QTLs 
detected were drawn using the MapChart 2.2 software 
[81]. 

Abbreviations 

AFLP: Amplified fragment length polymorphic; EC: Expansion coefficient; 
E-QE: E-QTL x Environment effect; E-QTL: Epistatic QTL; LD: Long-day; 
LGs: Linkage groups; LOD: Logarithm of odds ratio; MAS: Marker assisted 
selection; MCIM: Mixed-model based composite interval mapping; 
MQM: Multiple QTL model mapping; PDI: Popping dimension index; 
PUS: Percentage of unpopped seeds; QE: QTL x Environment effect; 
QTL: Quantitative trait loci; RFLP: Restriction fragment length polymorphic; 
RILs: Recombinant inbred lines; SD: Short-day; SNP: Single nucleotide 
polymorphism; SSR: Simple sequence repeats. 

Authors' contributions 

FY-L carried out genetic and molecular marker analyses, developed the 
genetic map and drafted the manuscript. MS developed the recombinant 
inbred line population, performed the field trials and carried out the 
phenotypic evaluation of plant material, and collaborating in writing the 
manuscript. CC and MG-A collaborated in marker genetic analyses and 
contributed to data analysis. MDLF collaborated in the phenotypic analysis of 
segregating population and the selection of microsatellite markers. JC 
participated in the design of genotyping experiments, supported mapping 
methodologies and contributed to a critical review of the manuscript. AMDR 
provided significant information about the plant material, collaborated in the 
analysis of the phenotypic and genetic data and reviewed the manuscript. 
RL conceived the project, planned the research work, assisted in analysis and 
interpretation of the data, and edited the manuscript. All authors have read 
and approved the final manuscript. 

Acknowledgements 

The authors thank Quival-Frutos Secos El Nogal (Pontevedra, Spain) for 
technical support and Diputacion de Pontevedra for farm facilities. We also 
thank Rosana Pereira Vianello Brondani from Embrapa Arroz e Feijao, CNPq 
(Brasil) for supplying some microsatellite primers. MDLF was supported by a 



research contract of the Xunta de Galicia. This work has been funded by 
grants PET2008_0167, EUI2009-04052 and AGL201 1-25562 of the Ministerio 
de Ciencia e Innovacion and PGIDI03RAG16E of the Xunta de Galicia. 

Author details 

^epartamento de Biologia Aplicada (Genetica), Edificio CITE ll-B, Centro de 
Investigacion en Biotecnologia Agroalimentaria (BITAL). Campus de 
Excelencia Internacional Agroalimentario, CeiA3, Universidad de Almena, 
Carretera de Sacramento s/n, 04120 Almena, Spain. 2 Departamento de 
Recursos Fitogeneticos, Mision Biologica de Galicia-CSIC, P.O. Box 28, 36080 
Pontevedra, Spain. 

Received: 14 March 2012 Accepted: 3 August 2012 
Published: 9 August 2012 



References 

1. Tohme J, Toro O, Vargas J, Debouck DG: Variability in Andean Nuna 
common beans {Phaseolus vulgaris, Fabaceae). Econ Bot 1995, 49:78-95. 

2. Gepts P, Osborn T, Rashka K, Bliss F: Phaseolin-protein variability in wild 
forms and landraces of the common bean {Phaseolus vulgaris L): 
evidence for multiple centers of domestication. Econ Bot 1986, 
40:451-468. 

3. Voysest O: Voriedodes de frijol en America Latino y su origen. Cali, Colombia: 
Spanish Centro Internacional de Agricultura Tropical; 1983. 

4. Kornegay J, White JW, Dominguez JR, Tejada G, Cajiao C: Inheritance of 
photoperiod response in Andean and Mesoamerican common bean. 

Crop So 1993, 33:977-984. 

5. CIAT: Morfologfa de la planta de frijol comun (Phaseolus vulgaris L). Guia de 
estudio. Cali, Colombia: Spanish. Centro Internacional de Agricultura 
Tropical; 1984. 

6. Hoseney RC, Zeleznak K, Abdelrahman A: Mechanism of popcorn popping. 
J Cereal Sci 1983, 1:43-52. 

7. da Silva WJ, Vidal BC, Martins MEQ, Vargas H, Pereira C, Zerbetto M, Miranda 
LCM: What makes popcorn pop. Nature 1993, 362:417. 

8. Li YL, Dong YB, Niu SZ, Cui DQ: QTL for popping characteristics in 
popcorn. Plant Breed 2007, 126:509-514. 

9. Spaeth SC, Debouck DG, Tohme J, van Beem J: Microstructure of nunas: 
Andean popping beans {Phaseolus vulgaris L). Food Microstruct 1989, 
8:263-269. 

10. Vorwald J, Nienhuis J: Effects of seed moisture content, cooking time, and 
chamber temperature on nuha bean {Phaseolus vulgaris L.) popping. 

HortScience 2009, 44:135-137. 

11. van Beem J, Kornegay J, Lareo L: Nutritive value of the nuna popping 
bean. Econ Bot 1992, 46:164-170. 

12. Kmiecik K, Nienhuis J: Development of nuna beans with temperate zone 
adaptation. Annu Rpt Bean Improv Coop 1 997, 40:36-37. 

13. Ogg B, Brick MA, Pearson C: Introgression of popping ability into dry 
beans adapted to Colorado. Annu Rpt Bean Improv Coop 1998, 41:35-36. 

14. Dofing S, D'Croz-Mason MN, Thomas-Compton MA: Inheritance of 
expansion volume and yield in two dent corn x popcorn crosses. Crop 
Sci 1991, 31:715-718. 

15. Ashman RB: Popcorn. In Cooperative extension service, Plant disease control, 
Bulletin BP-4. West Lafayette: Purdue University; 1984. 

16. Clary G: A study of the inheritance of expansion in popcorn. West Lafayette, 
IN: PhD thesis. Purdue University; 1954. 

17. Zeigler KE: Popcorn. In Specialty Corns. 2nd edition. Edited by Hallauer AR. 
Boca Raton, FL, USA: CRC Press; 2001:199-234. 

18. Li YL, Wu SW, Niu SZ, Li ZQ: Analysis of combining ability of popcorn with 
normal corn inbreds belonging to different heterotic groups. J Henan 
Agric Univ 2003, 37:1-5. In Chinese. 

19. Dofing SM, Thomas-Compton MA, Buck JS: Genotype x popping method 
interaction for expansion volume in popcorn. Crop Sci 1990, 30(1 ):62-65. 

20. Lu HJ, Bernardo R, Ohm HW: Mapping QTL for popping expansion 
volume in popcorn with simple sequence repeat markers. Theor Appl 
Genet 2003, 106:423-427. 

21. Babu R, Nair SK, Kumar A, Rao HS, Verma P, Gahalain A, Singh IS, Gupta HS: 
Mapping QTLs for popping ability in a popcorn x flint corn cross. Theor 
Appl Genet 2006, 1 12:1392-1399. 

22. Dudley JW: Molecular markers in plant improvement: manipulation of 
genes affecting quantitative traits. Crop Sci 1993, 33:660-668. 



Yuste-Lisbona et al. BMC Plant Biology 2012, 12:136 
http://www.biomedcentral.eom/1 471-2229/1 2/1 36 



23. Tanksley SD, Nelson JC: Advanced backcross QTL analysis: a method for the 
simultaneous discovery and transfer of valuable QTLs from unadapted 
germplasm into elite breeding lines. TheorAppI Genet 1996, 92:191-203. 

24. Bennett MD, Leitch IJ: Nuclear DNA amounts in angiosperms. Ann Bot 

1995, 76:113-176. 

25. Freyre R, Skroch PW, Geffroy V, Adam-Blondon AF, Shirmohamadali A, 
Johnson W, Llaca V, Nodari R, Pereira P, Tsai SM, Tohme J, Dron M, Nienhuis 
J, Vallejos C, Gepts P: Towards an integrated linkage map of common 
bean. 4. Development of a core linkage map and alignment of RFLP. 
TheorAppI Genet 1998, 97:847-856. 

26. Hanai LR, Santini L, Aranha LEC, Pelegrinelli MHF, Gepts P, Tsai SM, Carneiro 
ML: Extension of the core map of common bean with EST-SSR, RGA, 
AFLP, and putative functional markers. Mol Breeding 2010, 25:25-45. 

27. Broughton WJ, Hernandez G, Blair M, Beebe S, Gepts P, Vanderleyden J: 
Beans {Phaseolus spp.) - model food legumes. Plant Soil 2003, 252:55-128. 

28. Kelly JD, Gepts P, Miklas PN, Coyne DP: Tagging and mapping of genes 
and QTL and molecular marker-assisted selection for traits of economic 
importance in bean and cowpea. Field Crops Res 2003, 82:135-154. 

29. Kolkman JM, Kelly JD: QTL conferring resistance and avoidance to white 
mold in common bean. Crop Sci 2003, 43:539-548. 

30. Frei A, Blair MW, Cardona C, Beebe SE, Gu H, Dorn S: QTL mapping of 
resistance to Thrips polmi Karny in common bean. Crop Sci 2005, 
45:379-387. 

31 . Cichy KA, Blair MW, Galeno C, Snapp SS, Kelly JD: QTL analysis of root 
architecture traits and low phosphorus tolerance in an Andean bean 
population. Crop Sci 2009, 49:59-68. 

32. Blair MW, Medina Jl, Astudillo C, Rengifo J, Beebe SE, Machado G, Graham R: 
QTL for seed iron and zinc concentration and content in a 
Mesoamerican common bean {Phaseolus vulgaris L.) population. Theor 
Appl Genet 2010, 121:1059-1070. 

33. Blair MW, Astudillo C, Rengifo J, Beebe SE, Graham R: QTL analyses for seed 
iron and zinc concentrations in an intra-genepool population of Andean 
common beans {Phaseolus vulgaris L). Theor Appl Genet 201 1, 
122:511-521. 

34. Blair MW, Galeano CH, Tovar E, Muhoz-Torres MC, Castrillon AV, Beebe SE, 
Rao IM: Development of a Mesoamerican intra-genepool genetic map for 
quantitative trait loci detection in a drought tolerant x susceptible 
common bean {Phaseolus vulgaris L.) cross. Mol Breeding 2012, 29:71-88. 

35. Galeano CH, Fernandez AC, Franco-Herrera N, Cichy KA, McClean PE, 
Vanderleyden J, Blair MW: Saturation of an intra-gene pool linkage map: 
towards a unified consensus linkage map for fine mapping and synteny 
analysis in common bean. PLoS One 2011, 6:12. 

36. Welsh W, Bushuk W, Roca W, Singh SP: Characterization of agronomic 
traits and markers of recombinant inbred lines from intra- and interracial 
populations of Phaseolus vulgaris L. Theor Appl Genet 1995, 91:169-177. 

37. Gonzalez AM, Rodiho AP, Santalla M, De Ron AM: Genetics of intra-gene 
pool and inter-gene pool hybridization for seed traits in common bean 
{Phaseolus vulgaris L.) germplasm from Europe. Field Crops Res 2009, 
112:66-76. 

38. Koinange EMK, Gepts P: Hybrid weakness in wild Phaseolus vulgaris L 
JHered]992, 83:135-139. 

39. Singh SP, Gutierrez JA: Geographical distribution of the DL1 and DL2 
genes causing hybrid dwarfism in Phaseolus vulgaris L, their association 
with seed size, and their significance to breeding. Euphytica 1984, 
33:337-345. 

40. Xing YZ, Tan YF, Hua JP, Sun XL, Xu CG, Zhang QF: Characterization of the 
main effects, epistatic effects and their environmental interactions of 
QTLs on the genetic basis of yield traits in rice. TheorAppI Genet 2002, 

105:248-257. 

41. Kulwal PL, Singh R, Balyan HS, Gupta PK: Genetic basis of pre-harvest 
sprouting tolerance using single-locus and two-locus QTL analyses in 
bread wheat. Funct Integr Genom 2004, 4:94-101. 

42. Shen X, Zhang T, Guo W, Zhu X, Zhang X: Mapping fiber and yield QTLs 
with main, epistatic and QTL x environment interaction effects in 
recombinant inbred lines of upland cotton. Crop Sci 2006, 46:61-66. 

43. Yang X, Guo Y, Yan J, Zhang J, Song T, Rocheford T, Li JS: Major and minor 
QTL and epistasis contribute to fatty acid compositions and oil 
concentration in high-oil maize. TheorAppI Genet 2009, 120:665-678. 

44. Pedrosa-Harand A, Porch T, Gepts T: Standard nomenclature for common 
bean chromosomes and linkage groups. Annu Rep Bean Improv Coop 
2008,51:106-107. 



Page 15 of 16 



45. Blair MW, Buendia HF, Giraldo MC, Metais I, Peltier D: Characterization of 
AT-rich microsatellites in common bean {Phaseolus vulgaris L.). TheorAppI 
Genet 2008, 118:91-103. 

46. Blair MW, Hurtado N, Chavarro CM, Muhoz-Torres M, Giraldo MC, Pedraza F, 
Tomkins J, Wing R: Gene-based SSR markers for common bean {Phaseolus 
vulgaris L.) derived from root and leaf tissue ESTs: an integration of the 
BMc series. BMC Plant Bio 201 1, 1 1:50. 

47. Cordoba JM, Chavarro C, Schlueter JA, Jackson SA, Blair MW: Integration of 
physical and genetic maps of common bean through BAC-derived 
microsatellite markers. BMC Genomics 2010, 11:436. 

48. Perez-Vega E, Paheda A, Rodnguez-Suarez C, Campa A, Giraldez R, Ferreira 
JJ: Mapping of QTLs for morpho-agronomic and seed quality traits in a 
RIL population of common bean {Phaseolus vulgaris L). Theor Appl Genet 
2010, 120:1367-1380. 

49. Campos T, Oblessuc PR, Sforca DA, Cardoso JMK, Baroni RM, Sousa ACB, 
Carbonell SAM, Chioratto AF, Garcia AAF, Rubiano LB, Souza AP: Inheritance 
of growth habit detected by genetic linkage analysis using 
microsatellites in the common bean {Phaseolus vulgaris L.). Mol Breeding 
2011,27:549-560. 

50. Garcia RAV, Rangel PN, Brondani C, Martins WS, Melo LC, Carneiro MS, Borba 
TCO, Brondani RPV: The characterization of a new set of EST-derived 
simple sequence repeat (SSR) markers as a resource for the genetic 
analysis of Phaseolus vulgaris. BMC Plant Bio 201 1 , 1 2:41 . 

51 . Blair MW, Diaz JM, Hidalgo R, Diaz LM, Duque MC: Microsatellite 
characterization of Andean races of common bean 
{Phaseolus vulgaris L). TheorAppI Genet 2007, 1 16:29-43. 

52. Blair MW, Muhoz-Torres M, Giraldo MC, Pedraza F: Development and 
diversity assessment of Andean-derived, gene-based microsatellites for 
common bean {Phaseolus vulgaris L). BMC Plant Bio 2009, 9:100. 

53. McElroy JP, Zhang W, Koehler KJ, Lamont SJ, Dekkers JC: Comparison of 
methods for analysis of selective genotyping survival data. Genet Sel Evol 
2006, 38:637-655. 

54. Mackay I, Powell W: Methods for linkage disequilibrium mapping in 
crops. Trends Plant Sci 2007, 12:57-63. 

55. Ravi K, Vadez V, Isobe S, Mir RR, Guo Y, Nigam SN, Gowda MVC, 
Radhakrishnan T, Bertioli DJ, Knapp SJ, Varshney RK: Identification of 
several small main-effect QTLs and a large number of epistatic QTLs for 
drought tolerance related traits in groundnut {Arachis hypogaea L). 
Theor Appl Genet 201 1 , 1 22:1 1 1 9-1 1 32. 

56. Vorwald J, Nienhuis J: Heritability and correlations among popping 
characteristics and seed weight in a broad-based temperate-adapted 
nuna bean {Phaseolus Vulgaris L.) population. Hortscience 2009, 
44:214-216. 

57. Gu W, Zhu J, Wallace DH, Singh SP, Weeden NF: Analysis of genes 
controlling photoperiod sensitivity in common bean using DNA markers. 

Euphytica 1998, 102:125-132. 

58. SAS Institute Inc: SAS version 9.02. Cary, NC, USA: SAS Institute Inc; 2010. 

59. Doyle JJ, Doyle JL: Isolation of plant DNA from fresh tissue. Focus 1990, 
12:13-15. 

60. Vos P, Hogers R, Bleeker M, Reijans M, Van de Lee T, Hormes M, Frijters A, 
Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA- 
fingerprinting. Nucleic Acids Res 1995, 23:4407-4414. 

61. Noguera FJ, Capel J, Alvarez Jl, Lozano R: Development and mapping of a 
codominant SCAR marker linked to the andromonoecious gene of 
melon. TheorAppI Genet 2005, 1 10:714-720. 

62. Benchimol LL, Campos T, Carbonell SAM, Colombo CA, Chioratto AF, 
Formighieri EF, Gouvea LRL, Souza AP: Structure of genetic diversity 
among common bean {Phaseolus vulgaris L.) varieties of Mesoamerican 
and Andean origins using new developed microsatellite markers. Genet 
Resour Crop Evol 2007, 54:1 747-1 762. 

63. Cardoso JMK, Oblessuc PR, Campos T, Sforc 5 a DA, Carbonell SAM, Chioratto 
AF, Formighieri EE, Souza AP, Benchimol LL: New microsatellite markers 
developed from an enriched microsatellite common bean library. Pesq 
Agropec Bras 2008, 43:929-936. 

64. Gaitan-Solis E, Duque MC, Edwards KJ, Tohme J: Microsatellite repeats in 
common bean {Phaseolus vulgaris): isolation, characterization, and cross- 
species amplification in Phaseolus ssp. Crop Sci 2002, 42:2128-2136. 

65. Blair MW, Muhoz-Torres M, Pedraza F, Giraldo MC, Buendia HF, Hurtado N: 
Development of microsatellite markers for common bean {Phaseolus 
vulgaris L.) based on screening of non-enriched small insert genomic 
libraries. Genome 2009, 52:772-787. 



Yuste-Lisbona et al. BMC Plant Biology 201 2, 1 2:1 36 Page 1 6 of 1 6 

http://www.biomedcentral.eom/1 471-2229/1 2/1 36 



66. Blair MW, Pedraza F, Buendia HF, Gaitan-Solis E, Beebe SE, Gepts P, Tohme J: 
Development of a genome-wide anchored microsatellite map for 
common bean {Phaseolus vulgaris L). Theor Appl Genet 2003, 
107:1362-1374. 

67. Yu K, Park SJ, Poysa V, Gepts P: Integration of simple sequence repeat 
(SSR) markers into a molecular linkage map of common bean {Phaseolus 
vulgaris L). J Hered 2000, 91 :429-434. 

68. Buso GSC, Amaral ZPS, Brondani RPV, Ferreira ME: Microsatellite markers 
for the common bean Phaseolus vulgaris. Mol Ecol Notes 2006, 6:252-254. 

69. Grisi MCM, Blair MW, Gepts P, Brondani C, Pereira PAA, Brondani RPV: 
Genetic mapping of a new set of microsatellite markers in a reference 
common bean {Phaseolus vulgaris) population BAT93 x Jalo EEP558. 
Genet Mol Res 2007, 6:691-706. 

70. Hanai LR, Campos T, Camargo LEA, Benchimol LL, Souza AP, Melotto M, 
Carbonell SAM, Chioratto AF, Consoli L, Formighieri EF, Siqueira MVBM, Tsai 
SM, Vieira MLC: Development, characterization, and comparative analysis 
of polymorphism at common bean SSR loci isolated from genie and 
genomic sources. Genome 2007, 50:266-277. 

71. Galeano CH, Fernandez AC, Gomez M, Blair MW: Single strand 
conformation polymorphism based SNP and Indel markers for genetic 
mapping and synteny analysis of common bean {Phaseolus vulgaris L). 
BMC Genomics 2009, 1 0:629. 

72. Hougaard BK, Madsen LH, Sandal N, Moretzsohm MC, Fredslund J, Schauser 
L, Nielsen AM, Rohde T, Sato S, Tabata S, Bertioli DJ, Stougaard J: Legume 
anchor markers link syntenic regions between Phaseolus vulgaris, Lotus 
japonicus, Medicago truncatula and Arachis. Genetics 2008, 179:2299-2312. 

73. McConnell M, Mamidi S, Lee R, Chikara S, Rossi M, Papa R, McClean P: 
Syntenic relationships among legumes revealed using a gene-based 
genetic linkage map of common bean {Phaseolus vulgaris L). Theor Appl 
Genet 2010, 121:1103-1116. 

74. Montgomery J, Wittwer CT, Palais R, Luming Z: Simultaneous mutation 
scanning and genotyping by high-resolution DNA melting analysis. 
Nature 2007, 2:59-66. 

75. van Ooijen JH: JoinMap®4, software for the calculation of genetic linkage 
maps in experimental populations. Wageningen, Netherlands: 

Kyazma BV; 2006. 

76. Kosambi DD: The estimation of map distances from recombination 
values. Ann Eugen 1944, 12:172-175. 

77. van Ooijen JH: MapQTL®5, software for the mapping of quantitative trait loci 
in experimental populations. Wageningen, Netherlands: Kyazma BV; 2004. 

78. Yang J, Hu CC, Hu H, Yu RD, Xia Z, Ye XZ, Zhu J: QTLNetwork: mapping 
and visualizing genetic architecture of complex traits in experimental 
populations. Bioinformatics 2008, 24:721-723. 

79. Wang CS, Rutledge JJ, Gianola D: Bayesian analysis of mixed linear models 
via Gibbs sampling with an application to litter size in Iberian pigs. Genet 
SelEvol 1994, 26:91-115. 

80. Miklas PN, Porch T: Guidelines for common bean QTL nomenclature. 
Annu Rep Bean Improv Coop 2010, 53:202-204. 

81 . Voorrips RE: MapChart: software for the graphical presentation of linkage 
maps and QTLs. J Hered 2002, 93:77-78. 



doi:1 0.1 186/1471-2229-12-136 

Cite this article as: Yuste-Lisbona et al:. Marker-based linkage map of 
Andean common bean {Phaseolus vulgaris L.) and mapping of QTLs 
underlying popping ability traits. BMC Plant Biology 2012 12:136. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



